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LauremModulatoPPulses.nb 



" Needs [ "Laurent Funct ions * "] 

RuleDelayed: : rhs : Pattern t_ appears on the right-hand side of rule 

PhaseAngle[L_] [t_] (PhaseAngle[L] [t_] = Module[{xl, x2, x3 ( x4, x5, x6}, «1»J). 

Needs [ n LaurentNotationTest * " ] 

Needs: :nocont : Context LaurentNotationTest * was not created when Needs was evaluated. 

We develop the modulating function that is used in GSM first to use as an example to demonstrate Laurent's idea. 
3 

T := 

812500 
BT : = 0.3 

1 

Modulationlndex : = — 
2 

See LaurentTheory.nb to see the derivation of the C functions 

When running the notebook we fix the value of L here and the rest of the graphs and functions use this value 

" jfi'^i' — - <t 



■ % 4 , iu:r &-< 
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L := 6; 

PhaseFunction[t„] = N[$ x ,, t ]; 

phasepoints = Table [{t, PhaseFunction [t T] } // N, {t, 0, L , 1/40}]; 
ListPlot [phasepoints, Plot Joined -> True] 

NIntegrate: : ncvb : 

NIntegrate failed to converge to prescribed accuracy after 7 recursive bisections 
in LaurentFunctions * Private * tl near LaurentFunct ions * Private % tl = -1 . 91796 x 10" 7 . 

NIntegrate: :ncvb : 

NIntegrate failed to converge to prescribed accuracy after 7 recursive bisections 
in LaurentFunctions * Private* tl near LaurentFunct ions * Private * tl = -5 . 52869 x 10 -7 . 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method->Oscillatory in NIntegrate. 

NIntegrate: :ncvb : 

NIntegrate failed to converge to prescribed accuracy after 7 recursive bisections 
in LaurentFunctions* Private* tl near LaurentFunctions * Private * tl = -1 . 96696 x 10" 7 . 

General :: stop : 

Further output of NIntegrate: : ncvb will be suppressed during this calculation. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method- >Oscillatory in NIntegrate. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method- >Oscillatory in NIntegrate. 

General : : stop : 

Further output of NIntegrate: : slwcon will be suppressed during this calculation. 




- Graphics 



BitSeq = Table [1, {i, 1, 20}] 

{1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, i f i, l f 1} 

RandomBitSeq = Table [Random [Integer, {0, 1}], {i, 1, 40}] //Map[# (-2) + It, #] & 

{i, i, -l, -l, -i, i, i, -i, i, i, i, i, i, -i, i, -i, i, i, i/ _ 1/ _ lt _ x 
-1, 1, 1, -1, -1, -1, -1, 1, 1, -1, -1, 1, -1, -1, -1, -1, 1} 

Testing The Modulator ... 
Calculating the bandwidth of the signal 

xll = Table [LaurentC [L] [0] [tT], {t, 0, L + 1 - 1/200, 1/200}]; 

((Drop[xll, -1] -Drop[xll, 1]) 200 ) A 2 / 200 // Apply [Plus, #]& 

1.29684 



Here the pulse 6T is just stretched to 8T to see the performance 
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3 



'xample of stretching a C 0 pulse Of 6 to 8 

In this section we consider filtering the pulse 

fil = BesselFilter [8] [LowPass3dB [ 1 / T 2 Pi 0 . 3] ] [FrequencyResponse] [s] ; 

Plot [20 Log[10, Evaluate [ (fil /. s -> 2 Pi I f / T) //Abs]], {f, 0, 2}, 
GridLines -> Automatic, AxesLabel -> {"1/T", "dB"}] 



dB 



-20 



-40 



-60 



-80 



1/T 



- Graphics 



Res[t_] = InverseLaplaceTransf onn[f il, s, t] 



2027025 ((9.36812 + 0. I) E" 897172 -* Cos[139301. t] - 

(13.7944 + 0. I) E" 835672 - t Cos[420044. t] + (4.73111 + 0. I) E" 701358 - c Cos[708768. t] - 



(0.304883 + 0. I) E- 



Cos[1.02016xl0 b t] + 



(29.1264 + 0. I) E" 897172 - t Sin[139301. t] - 
(9.89398 + 0. I) E" 835672 ■ fc Sin[420044. t] - (0.321205 + 0. I) E' 701358 - c Sin[708768. t] + 
(0.375154 + 0. I) e"* 455818 * fc Sin [1 . 02016 x 10 6 t] ) 

Plot [Res [t T] , {t, 0, 6}, PlotRange ->'A11] 




1 2 3 4 5 6 

- Graphics - 

filTaps = Table [Res [t T] , {t, 0, 6, 1/20}] //Chop; 

We estimate the group delay to be 1 .8T 
Length [ f i ITaps ] 



121 
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L : = 8 ; 

PhaseFunction[t_] sN[$x., t ]; 

phasepoints = Table [{t, PhaseFunction[t T] } // N, {t, 0, L, 1/40}]; 
ListPlot [pbasepoints, Plot Joined -> True] 

NIntegrate: :ncvb : 

NIntegrate failed to converge to prescribed accuracy after 7 recursive bisections 
in LaurentFunctions * Private 4 tl near LaurentFunctions % Private * tl = - 3 . 68196 x 10" 7 . 

NIntegrate: :ncvb : 

NIntegrate failed to converge to prescribed accuracy after 7 recursive bisections 
in LaurentFunctions * Private 4 tl near LaurentFunctions 4 Private * tl = -7 . 4914 x 10~ 7 . 

NIntegrate: :ncvb : 

NIntegrate failed to converge to prescribed accuracy after 7 recursive bisections 
in LaurentFunctions * Private * tl near LaurentFunctions ' Private ' tl = - 5 . 43232 x 10" 7 . 

General: : stop : 

Further output of NIntegrate: :ncvb will be suppressed during this calculation. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method- >Oscillatory- in NIntegrate. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method->Oscillatory in NIntegrate. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method->Oscillatory- in NIntegrate. 

General :: stop : 

Further output of NIntegrate: : slwcon will be suppressed during this calculation. 
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Graphics - 



pulse = Table [ Laurent C [8] [0] [tT] ; {t, 0, 9, 1/20}]; 



pulse2 = Table [LaurentC [8] [1] [t T] , {t, 0, 7, 1/20}]; 
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FilteredPulse = T/20 MapConvolve [pulse, filTaps] //Take[#, {34, I,ength[#] }]& 



-15 
-13 
-11 



, 5.35887 xl0" 15 , 
, 9.65729xlCT 13 , 
, 8.83867 X10" 11 , 
4.38922 xlCT 9 , 8 
2.32319xl0" 7 , 4 
3 .75853 x 1CT 6 , 6.23xl0" 6 , 0 
0.000061746, 0.000093496, 0 



1.60413 x 10" 14 
2.51044 xl0~ 12 
2.02275x 10 
.95697 xl0" 9 
.20166 x 10" 7 



10 



4.66278 xl0" 14 
6.35507xl0" 12 
4.51958 xl0" 10 
1.78914 xl0" 8 , 3 
7.46062 xl0" 7 , 1 
0000101618, 
00013964, 0 



1.31688 x 10"", 
1. 56747 x 10" 11 , 
9.86455 xl0" 10 , 
4999xl0~ 8 , 6 .70842 x10" 
30124 x 10" 6 , 

0.0000163174, 0.0000258049, 
000205786, 0.000299337, 



{1.73731x10 
3.61578 x 10 
3.76889 x 10 
2.10426 xl0" 9 , 
1.26054 x 10" 7 , 
2 .23035 x 10' 6 , 
0.0000402068, 

0.000429924, 0.000609^5, 0.000854853, 0.00118424; 0.00162196, 0 00219697 
0.0029439, 0.00390362, 0.00512368, 0.00665873, 0.00857068, 0.0109288, 0 0138093 
.0172954, 0.021476, 0.0264451, 0.0323001, 0.0391407, 0.0470666, 0.0561754 
.0665603, 0.0783079, 0.0914948, 0.106186, 0.12243, 0.140261, 0.159691 0 180711 
.20329, 0.22737, 0.25287, 0.279682, 0.307673, 0.336686, 0.366539, 0.397029 
.427934, 0.459013, 0.490013, 0.520668, 0.550708, 0.579856, 0.607837, 0 634379 
.65922, 0.682107, 0.702805, 0.721098, 0.73679, 0.749713, 0.759727, 0 76672 
.770616, 0.771369, 0.768969, 0.763442, 0.754846, 0.743274, 0.728852, 0 711736 
0.670178, 0.646176, 0.62035, 0.592964, 0.56429, 0.534607, 0.504193 ' 
0.442277, 0.411305, 0.380656, 0.350558, 0.321221, 0.29283, 0 26555* 
.214852, 0.191636, 0.169936, 0.149791, 0.13122, 0.114217, 0.0987613 
.0612016, 0.0513967, 0.0428146, 0.035365, 0.0289542, 



0.692108, 
0.473326, 
0.239519, 0. 

0.0848118, 0.0723143, 0. 

0.023487, 0.0188687, 0.0150065, 0.0118109, 0.00919671, 0.00708412, 0 0053992 
0.00407439, 0.00304877, 0.00226811, 0.0016848, 0.00125761, 0.000951417 0 000736739 
0.000589309, 0.000489543, 0.000422011, 0.000374899, 0.000339478, 0.000309*6 
0.000281223, 0.000251974, 0.000220764, 0.000187441, 0.000152497, 0 000116824 
0.0000815146, 0.0000477044, 0.0000164566, -0.0000113198, -0.0000349206 
-0.0000538729, -0.0000679391, -0.0000771062, -0.0000815622, -0.0000816 651 
-0.0000779073, -0.00007 0877, -0.000061222, -0.0000496142, -0.0000367186, ' 
-0.0000231668, - 9 . 53584 x 10" 6 , 3 . 66817 x 10" 6 , 0.0000160205, 0.0000271838 
0.0000369093, 0.0000450345, 0.0000514769, 0.0000562268, 0.0000593363 
0.0000610882, 0.0000600456, 0.0000579702, 0.0000550589, 0 0000515076' 
0.0000432203, 0.000038813, 0.00003 44153, 0.0000301387, 0.0000260716 ' 



0.000060909, 
0.0000475036, 



0. 0000222801, 
8 .4469 x 10" 6 , 
2.3199xl0 -6 , 
4.45602 xl0" 7 , 
5.71866 xl0" 8 , 
4.63953 xl0" 9 , 

.2234 xl0" 10 , 
.92683 x 10' 12 , 
.95416xl0" 14 , 
.65283 xl0" 16 , 
-1.01228 xl0" 17 
1.9917xl0" 19 , 

1. 94753 xlO" 21 
7.07296 xlO" 23 



0.0000188096 
6. 69967 x 10~ 6 , 
1.71859xl0" 6 , 
3 .05722 x 10" 7 , 
3 .59686 xl0" 8 , 
2.64191xl0 -9 , 
1.12886 X 10" 10 , 
, 2.65259 xl0" 12 
, 3.14517xl0" 14 
, 1 . 93314 x 10" 16 
6.49877x10 
-1.14008* 10" 19 
3.72298xl0" 22 
,6.84862 xl0" 23 



•5.51867xl0" 25 , 7.52915 x 10' 



0.000015687, 
5.24481xl0" 6 , 
1.25466xlO" 6 , 
. 2.06318 xl0" 7 
2 .22017 x 10" 8 
1.4722xl0" 9 , 
, 5.59276X10" 11 , 
1.14156X10" 12 , 
1.23089xl0 r14 , 
9.87738 xlO" 17 , 
B , -2.08455 x 10 
-1.66045xl0" 20 
2.02535 xl0" 22 , 
, 3.68916xl0" 23 , 
26 , 5.81584 xl0 -25 



0.0000129228, 0.000010514, 
4.05145x 10" 6 , 3.08723 xl0" 6 



9.0236 xl0* 7 , i 
1.36898xlO" 7 
1.34418 x 10" 8 , 7 
8.02333 xl0" 10 , 4 
2.70717 X 10" 11 , 
4.73562 xl0" 13 
4.73302 x 10 
2 .31472 x 10 
18 , 1.80177x10 
, 3 .7877 x 10" 21 , 
1.36007 xl0" 22 , 
2 . 6013 x 10 



-8 

-9 
-10 



-15 
-18 



^-24 



39109 x 10" 
8 . 92721 x 10" 
97812 x10" 
27388x 10" 
1.28211x 10" 11 , 
. 1. 94329 x 10" 13 , 
1.66595 xl0" 15 , 
-1.50823 xlO" 17 , 
19 , 1.70062 x 10" 19 , 
1 .47752 x 10" 21 , 
6.72514 xl0" 24 , 
1.51251xl0" 24 , 



1.75371 x 10" 25 , 1.22392 xl0" 26 , 0} 



LaurentkodulatorPulses.nb 
* — 



l 6 



FilteredPulse2 = T/20 MapConvolve [pulse2, filTaps] //Take[#, {50, Length[#j )]& 

{1.27972X1CT 7 , 2.18306xl0" 7 , 3 . 66972 x 10" 7 , 6.0806xl0" 7 , 9 . 93427 x 1 CT 7 , 1.60076xlCT 6 
2.54471xl'0" 6 , 3.99199xl0" 6 , 6 . 18149 x 1CT 6 , 9 . 45059 x 10" 6 , 0.0000142688 0 0000212804 
n'^^^^f 2 '. 0 - 0000456572 ' 0*0000657079, 0.0000934812, 0 . 0001 14 93 0. So0182902 
S"22??^^ 6 ' 0.000342391, 0.000460925, 0.000613932, 0.000809185, 0 00105553 
0.00136282, 0.00174186, 0.00220418, 0.00276183, 0.00342705, 0 00421186 C l 00512761 
0.00618453, 0.00739101, 0.00875319, 0.0102744, 0.0119544 0 0137895 0 0157714 ' 
0.0178876, 0.0201213, 0.0224509, 0.024851, 0 . 0272 92 1 0 . 02 97419 0 03 2i652 
0.0345256, 0.0367857, 0.0389085, 0.0408581, 0.0426011, 0.044107 0*0453496 
0.0463073, 0.0469636ro;0473081, 0.0473359, 0.0470483, 0.0464525, 6 0455611 
°* 042967 ^ 0.0413131, 0.039459, 0.0374364, 0.0352785, 0.6330189,' 
0.0306913, 0.0283283, 0.0259611, 0.0236189, 0.021328, 0.019112, 0 0169912 
0.0149823, 0.0130987, 0.0113504, 0.0097438, 0.00828249, 0.00696697, 0 00579518 
0.00476282, 0.00386365, 0.0030899, 0.00243259, 0.00188191, 0.00142752 0 0010588S 
2-2222f55f ?4 °A°^" 6896 ' 0.000363597, 0.000236336, 0.000146695 "oOOOS^S ' 
S*222^S 7313 ' 0-0000318118, 0.0000253258, 0.0000271231, 0.000033 8422, 0 0000428466 
iiliS* 0.0000603451, 0.0000665053, 0.0000701193 0.0000710042 0 0000692325 
0.0000650647, 0.000058 88 85, 0.0000511665, 0.0000423925, 0.0000330554; 0 . 0000236122 ' 
0.0000144671, 5.95924xl0" 6 , - 1 . 64542 x 10' 6 , - 8 . 15591 x 10" 6 , -0.0000134543 
-0.0000174907, -0.00002 02751, -0.0000218687, -0.0000223734, -0.0000219217 
-0.000020666, -0.0000187698, -0.0000163988, -0.0000137138, -0.0000108643 ' 
-7.98457 X 10" 6 , - 5 . 18972 x 10" 6 , - 2 . 57416 x 10" 6 , - 2 . 1065 x 10* 7 , 1 . 84943 x 10' 6 ', 
3.57494xl0' 6 , 4 . 95319 x lO" 6 , 5 . 98756 x 10~ 6 , 6.69487 x10^, 7 . 10243 x 10" 6 , 
7.24522 xlO' 6 , 7 . 16306 x 10~ 6 , 6 . 89804 x 10" 6 , 6 . 4922 3 x 10" 6 , 5 . 98586 x 10* 6 , 
5.41573 xlO" 6 , 4.81423 xl0" 6 , 4 . 20863 x 10' 6 , 3.6208xl0" 6 , 3 . 0672 8 x 10* € , 2.55958xl0~ 6 
2.10471xl0" 6 , 1.70583 X 10' 6 , 1 . 3 62 97 x 10" 6 , 1 . 0737 6 x 10" 6 , 8 . 34162 x 1 O" 7 , 6.3908 xl0~ 7 ' 
4.82884xl0" 7 , 3 . 59 854 x 10" 7 , 2 . 64486 x 10" 7 , 1 . 91715 x 10" 7 , 1 . 37044 x 10~ 7 , 
9.65984 xl0' 8 , 6 . 71324 x 10' 8 , 4 . 59913 x 10" 8 , 3 . 10539 x 10~ 8 , 2 . 06612 x 10' 8 \ 1.3542 xl0" 8 
8.74155 xl0" 9 , 5.55593 xl0- 9 , 3 . 47564 x 10* 9 , 2 . 13956 x 10~ 9 , 1 . 29613 x 10~ 9 , 
7.73129X10" 10 , 4.54101xl0- 10 , 2 . 62 03 1 x lO' 10 , 1 . 4788 8 x 10" 10 , 8 . 15113 x 10" 11 , 
4.40884X10" 11 , 2.37341X10" 11 , 1 . 2 6865 x 10" 11 , 6 . 58958 x 10' 12 , 3 . 34992 x 10" 12 ,' 
1.6239 xl0* 12 , 7.01037 x 10" 13 , 2 . 33036 x 10" 13 , 7 . 75051 x 10" 14 , - 3 . 22249 x 10" 15 ', 
-8.12875xl0- 14 , -7.89371xl0" 14 , -5 . 29088 x 10" 14 , -3 . 14297 x 10~ 14 , - 8 . 21475 x 10" 15 
4.12557xl0" 15 , 1.9906 xl0- 15 , 1 . 82068 x 10' 15 , - 4. 3 6969 x 10" 15 , 4 . 67143 x 10" 16 , 
7.65162 xl0" 16 , 5.67707xl0" 16 , 4 . 93342 x 10" 16 , 7 . 52562 x lO' 17 , 7 . 75606 x 10' 17 
4.50401X10" 17 , -2.18469X 10" 17 , 1 . 39012 x 10" 16 , 1 . 92164 x lO' 16 , 6 . 31946 x 10~ 17 , 
-6.72328xl0' 18 , 1 . 247 04 x 10" 17 , 1 . 42063 x 10" 17 , 7 . 33335 x 10~ 18 , 1 . 13717 x 10' 19 
3.90929xl0 -19 , 1.30285xl0" 19 , 0} 

FiltPulse[8] [0] = {(Table[tT, {t, 0, 20, 1/20}] // N) // 

Take[#, Length [FilteredPulse] ]&, FilteredPulse} //Transpose// 
Interpolation 

InterpolatingFunctiont { {0, 0.0000492923}}, <>] 

FiltPulse[8] [1] = { (Table [t T, {t, 0, 20, 1/20}] // N) // 

Take[#, Length [FilteredPulse2] ]&, FilteredPulse2} //Transpose// 
Interpolation 

InterpolatingFunction[ { {0, 0.00003 89538}}, <>] 
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Plot[FiltPulse[8] [0] [tT], {t, 0, 9}, PlotRange -> All] 
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- Graphics - 

Plot [LaurentC [ 8] [0] [tT], {t, 0, 9}, PlotRange -> All]. 
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- Graphics - 
Show {%, %%] 
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Graphics - 



LaurentModulatorPulses.nb 



Plot [FiltPulse [8] [1] [t t] , {t, 0, 7}, PlotRange -> All] 
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- Graphics - 

Plot [ Laurent C [8] [1] [t T] , {t, 0, 7}, PlotRange -> All] . 
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- Graphics - 

Bandwidth [TestPulse_, Supp_] : = 
Module [{xll}, xll = Table [TestPulse [t T] , {t, 0, Supp - 1/200, 1/200}]; 
<(Drop[xll, -1] -Drop[xll, 1]) 200) *2 / 200 // Apply [Plus, #]&] 

BandWidth[FiltPulse[8] [0] , 9] 

0.711659 

Bandwidth [LaurentC[ 8] [0], 9] 

1 .29684 

Bandwidth [FiltPulse [8] [1] , 7] 

0.00338446 

torn = Modulator [L] [Table [1, {i, 1, 100}], SamplingXnterval -> T/10, 
NumberOf Curves -> 2, ModulatingPulse -* FiltPulse]; 



torn = Modulator [L] [RandomBitSeq, SamplingXnterval -* T/10, NumberOf Curves 
ModulatingPulse -> FiltPulse] ; 
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Li st Plot [Im [torn] , Plot Joined -> True] 




- Graphics - 

ListPlot [Re [torn] , Plot Joined -> True] 




D00 



- Graphics - 

ListPlot [Abs [torn] , Plot Joined -> True] 
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{Re[tom], Im[tom]} // Transpose // ListPlot [#, PlotJoined -> True, AspectRatio -> 1]& 




- Graphics - 

torn = Modulator [I,] [RandomBit Seq, Samplinglnterval -+ T/10, NumberOf Curves -> 1, 
ModulatingPulse -► FiltPulseJ ; 

Lis tPlot [Im[ torn] , PlotJoined -> True] 




000 



- Graphics - 
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ListPlot [Abs [torn] , Plot Joined -> True] 




- Graphics - 

{Re[toxn], lm[tom] } // Transpose // ListPlot [#> PlotJoined -> True, AspectRatio -> 1]& 




- Graphics - 

xll = Table[TestPulse[6][0][tT], {t, 0, 7 - 1/200, 1/ 200}]; 
((Drop[xll, -1] -Drop[xll, 1]) 200 )*2 / 200 // Apply [Plus, #]& 

0.92632 

TestPulse[6] [0] [t_] =LaurentC[4] [0] [t 5/7 ] 

1 ' Sl »['i 4 ' 8lAoF + J ^]] Si »W 4 < 40A5O ^]]sin[„[4 ( + 
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L := 5; 

PhaseFunction[t_] = #x, ft ; 

phasepoints = Table [{t, PhaseFunction[t T) } // N, {t, 0, L , 1/40}]; 
ListPlot [phasepoints, Plot Joined -> True] 

NIntegrate: : ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl - - 1 . 48053 x 10~ 6 . 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in-.-tl near tl - - 1 . 12255 x 10~ 6 . 

NIntegrate: : ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 08808 x 10" 6 . 

General : : stop : 

Further output of NIntegrate: : ncvb will be suppressed during this calculation. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration. is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method->Oscillatory in NIntegrate. 
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- Graphics - 
L := 6 

TestPulse2 [6] [0] [t_] = LaurentC [4] [0] [t 6/7] 



6 t 



Sin[#[4. JLL]] 



torn = Modulatory] [RandomBitSeq, Samplinglnterval T/40, NumberOf Curves -> 1, 
ModulatingPulse -* FiltPulse] ; 



Save [ "ModulatorDat a . m H , {T, L, RandomBitSeq, FiltPulse, torn}] 
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ListPlot [Xm[tom] , Plot Joined -> True] 
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- Graphics - 

ListPlot [Re [torn] , PlotJoined -> True] 
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- Graphics - 

ListPlot [Ai>s [torn] , Plot Joined -> True] 
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Graphics 
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{Re[tom], Im[tom]} // Transpose // Lis tPlot [#, PlotJoined -> True, AspectRatio -> l]& 




- Graphics - 

xll = Table[TestPulse2[6] [0] [t T] , {t, 0, 7 - 1/200, 1/200}]; 
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Needs ["Laurent Functions % "] 



RuleDelayed: :rhs : Pattern t_ appears on the right-hand side of rule 

PhaseAngle[L_] [t_] (PhaseAngle[L] [t_] = Module[{xl, x2, x3 ( x4 ( x5, x6}, «1» ) ) . 

Needs [ "LaurentNotationTest * n ] 



Needs: mocont : Context LaurentNotationTest' was not created when Needs was evaluated. 



Information on the functions used can be obtained using help. 



Name s [ " Laurent Func t i ons A * " ] 

{AKN, AlphaKI, ANKIni tialStateSetUp, BT, FiltPulse, h, hFiltered, Initialstate J 
LaurentC, LaurentLK; LaurentS, M, ModulatingPulse, Modulationlndex, Modulator 
NumberOf Curves, PhaseAngle, Receiver, Receiver Proper, S, Samplinglnterval 
StartingQuadrant, SyncSample, T, C, § , ij/} 

3 

X : = 

812500 
BT : = 0 - 3 



Modulationlndex : = — 
2 

<< ModulatorData.ro; 



RandomB i t S e q 



{1, 1, -1, -l f -1, 1, l, -l, l, -i # i, -i, i, _i, _x, .! . x ! x ! _ x , ± x 

1, -1, 1, -1, 1, -1, 1, 1, -1, -1, i, -l, 1, -i, x, 1, 1, -1, -1, -1, 1, -l 

i, i, i, -i, i, -i, -i, i, i, i # .i # .i, i, a, !, x, x, i, i, ± / lt : lf : x ,'i, _i, lt _ 1} 
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T 



3 



812500 



Snt+at 




^K t n'T + AT 



With L = 8, and M = 2 we can utilise only the main values. It so happens that when K = 1, the dominant value occur; 
2.5T in the filtered pulse. The three dominant values when K =0 occurs at 3.5T, 4.5T and 5.5T. The value of the pul: 
at 6.5T is smaller than the value of the pulse with K = 1 at 2.5 T. Thus the expression becomes 



S N rr«. iT - J**-*- 1 Pulse [0] [3T + <5T] + j^o.h-* Pulse [0] [4 T + <5T] + 
j*o,ii-5 p u ise [0] [5T + <5T] + J*!**-* Pulse [1] [2 T + <5T] 

Sht + at = J Ao ' B - , Pulse[0][3T + <5T] + J* 0 '"-* Pulse[0][4T + <5T] + 
j*o.»-s Pulse[0][5T + <5T] + j( a o-s-i^h-3) Pulse[l][2T + <5T] 

Given that a N - 3 has been already decided, we can precalculate the possibilities and store them. 



Sht*at= J* 0 -*-* ( j< a »-4 + «»-3> Pulse[0][3T + <5T] + J*«-* Pulse[0][4T + 6T] + 
Pulse [0] [5 T + 6T] + J< B »- * * °n- i ♦ «n- 3 * ««- > > Pulse [1] [2 T + 5T] ) 



S» t + at = J* 0 '"" 5 ( + Pulse[0][3T + <5T] + J**-< Pulse[0][4T + <5T] + 

Pulse [0] [5T + <5T] + Jt a »- * ♦ a «-a ) pulse [1] [2 T + 6T] ) 

IiOokUpTable[Pulse_, 6T_] : = 
Module [{xl, x2, x3, x4}, 

Xl = {{-1, -1, {-1, -1, 1}, 

{-1, 1, -1}, {-1, 1, 1}, {1, -1, -1}, {1, -1, 1}, {i/ 1, {1# X/ 

x2 [bitseq_] : = 

( J bit S e< I t[3]3*bit e a Q Et2]) Pulse[0][3T + 6T] + j*>^°*<*tl*U pui S e[03[4T + <5T] + 

Pulse[0][5T + <5T] + o toitoe<I ^ 3 n ♦bAfq[[i)j pulse[l][2T + <5T] ) ; 
{xl, Map[x2[#]&, xl]} // Transpose] 

tab = LookUpTable[FiltPulse[8] , 0.5 T] 

{{{-1, -1, -1}, -0.00122183 - 0.770616 1}, {{-1, -1, 1}, 0 . 70233 9 + 0 . 770 616 I } , 
{{-1, 1, -1}, 0.614125 - 0.770616 1}, {{-1, 1, 1}] 0.0869922 + 0.770616 1}, 
{{1, -1, -1}, 0.0869922 - 0.770616 1}, {{1, -1, 1}, 0 . 614125 + 0 .-770616 1} , 
{{1, 1, -1}, 0.702339 - 0.770616 1}, {{1, 1, 1}, -0.00122183 + 0.770616 1}} 

Sort[tab, #1[[1, 3]] >#2[[1, 3]]&] 

{{{1, 1, 1}, -0.00122183 + 0.770616 1}, {{1, -1, 1}, 0.614125 + 0.770616 1}, 
{{-1, 1, 1}, 0.0869922 + 0.770616 1}, {{-1, -1, 1}, 0.702339 + 0.770616 1}, 
{{1, 1, -1}, 0.702339 - 0.770616 1}, {{1, -1, -1}, 0.0869922 - 0.770616 1}, 
{{-1, 1, -1}, 0.614125 - 0.770616 1}, {{-1, -1, -1}, -0.00122183 - 0.770616 1}} 

FxltPulse [8] [0] [4 T + 0.5 T] 




But J 2 = 1 and so we get 



0.770616 



Now we build a receiver! ! First generate the modulated sequence 
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Options [Modulator] 

{startingQuadrant -> 0, 
InitialState-* {1, 1, 1, 1, l, i, l, i, i, i, i, i, i, lt lf x x lf x 1} 

3 

Samplinglnterval 2 6000000 ' Number ° f Curves -> 4 , ModulatingPulse Laureate} 



26000000 



1 

32 



tom2 = Modulator [L] [RandomBitSeq, NumberOf Curves -+ 2, ModulatingPulse FiltPulse] ; 
ListPlot [{Re[tom2] , Iro [tom2] } //Transpose, PlptJoined -> True] 

- Graphics - 
RandomBitSeq 

{1 ; 1 ; " x ' - 1 ' x ' 1 > x ' !< i< -i* -i, -1, -1, 1, 1, 1, -1, -1, -i, 1,11 
i, i, i # -i, 1, -1,1, i; 1, -i # _i, lf i # X/ lt lt lf lt _ lf _^ X/ lt _ 1} 

Receiver [L] [tom2, StartingQuadrant -> 4] 

{1, 1, 1, 1, 1, 1, -1, -1, -1, 1, l, -1, i, i, _ 1# lt _ x _! _ a _ x ii!.! , 

-1 1 1, 1, -1, 1, -1, -11, 1, 1, x, i, !, _ 1# 1# J 1§ 2 . x X ' 

x : 1 ' 1 1 ' 1 ' x ; 1- 1< -1. -1. 1. -1. 1, -1, 1. 1, 1, 1. -1, 

-1, -1, 1, -1, 1, 1, 1, -1, 1, -1, -1, 1, l; 1, _1, _ 1# 1, lf lf lf lt lf lt ' lf lf _ 1} 

We now calculate the BER given the |j- and L = 8. We write all the possible terms 

pulB.[0][T+«T] + J (a H .^a s . J+ a B . l)pulse[0][2T + i5T] + 

j(aB.4 + a H . 3 ) puisetOHST + 5 T ] + j*"-* Pulse[0][4T + <5T] + Pulse[0][5T + <5T] + 
jf-«N. 5 ) p u lse[0][6T + <5T] + j<- a »- Pulse[0][7T + <5T] .+ 

j(-«H-5 " a "-< ~ a H-7 ) Pulse [0] [8 T + <5T] + j(»N-4*«M-3 ♦ •» . 2 * « M - 1 ♦ «» - « H . 1 > Pul S6 [ 1 ] [ $T ] + 

j<.»-4*...->. -« H _ 2 > p u lse[l] [T + <5T] + J^.4^.|.^^.») PulseTl] [2 T + <5T1 + 

j<«»-4-" N - 3 ) Pulse[l][3T + <5T] + jCn-4-..-s > Pulse[l][4T + <5T] + 
j(-n-c) Pulse[l][5T + <5T] + j(-««.s^«.7 ) Pulse[l][6T + <5T] ) 

We select the imaginary 

0*..-. ( j^-*—^^-— pulse[0] [6T] + — 1 + « M . 3)pulfie[0][2T + 5T] + 

CT 3 "-* Pulse[0][4T + <5T] + ,T<- a »-5) Pul se [0][6T + 6T] + 

J (- N .5-a N . ( -a 1) . 7)pulse[0][8T + (5T] + .^.^ p ulse[ l] [T + 5T] + 

j(«n-4*« n -3 -«n-4 > p u lse[l][3T + <5T] + jt-«N. fi -> Pulse[l][5T + <5T] ) 

ModulationValue [Pulse_] [ {x0_, xl_, x2_, x3_ # x4_, x5„, x6_, x7_) ] [5T_] : = 
(1. * i)«o«i«a«*«4 ^[oj [<5T] + (1- , 1)xa+3 c3 + x4 pulse[0] [2 T + OT] + 

I)* 4 Pulse[0] [4T + 5T] + (1. * X)' x5 Pulse[0] [6T + <5T] + 
(1.* I)-**-*-' Pulse [0] [8T + 5T] + <1. » I)* 1 ****-* Pulse [1] [T + <5T] + 
X) xZ Pulse[l] [3T + 5T] + (1.- I)"* 6 Pulse [1] [5 T + <5T] 
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ModulationValueEFiltPulsetS]] [{-1, -1, -1, -l, -i, -i, _ lr -1>][0.5TJ // N 

-0.754091 I 

ConvertToBitSeq[Depth_] [Num_] : = 
Module [{bit, n} , 
n = Num; 

ConvertNext := (bit = Mod[n, 2];n = Floor [n / 2] ; bit) ; 
( Table [ConvertNext, {i, 1, Depth}] //Reverse) /. {0 -1}]? 

v •■• • - 

ConvertToBitSeq[8] [2] 

{-1, -1, -1, -1, -1, -1, 1, -1} 

Table [i, {i, 1, 2*8}] // Map [ConvertToBitSeq[8] , #]&; 

Select [Table [i, {i, 1, 2*8}] // Map [ConvertToBitSeq[8] , #]&, (#[[5]] == l)fi] 



{ 



-1, -1, -1, -1, 1, -1, -1, -1), 

-1, -1, -1, -1, 1, -1, -1, 1}, {-l, -l, _i, x, -1,-1, -1}, 

-1, -1, -1, -1, 1, -1, 1, 1}, {-1, -X, -1, -1, 1, 1, -l, -X}, 

-1, -1, -1, -1, 1, 1, -1, 1}, {-1, -X, -1, -X, 1, 1, 1, -1}, 

-1. -1, -1, -1, 1, 1, 1, 1}, {-1, -1, -X, 1, 1, -1, -1, -1}, 

-1, -1, -1, 1, 1, -1. -1, 1}, {-1, -X, -1, 1, 1, -1, 1, -1}, 

-1» -1' !. 1 - !. u. -1. -1' 1. 1. 1. -1. -1}. {-1. -1. -1, 1, 1, 1, -1, 1}, 

-1, -1, 1, -1, 1, -1, -1, 1}, {-1, -x, 1, -1, 1, -1, 1, _x}. 

"I' "J* 1- "I- !• 7 1 ' 1' J-)' _1 ' 1 ' _1 ' 1 ' X ' - 1 ' "I- !• "I' !• 1. "I' 1>. 

-}' }■ 1' 1 ' -!>' - 1 ' X ' - 1 ' X » 1 ' X >» t" 1 ' "I- 1. 1' I- "I. "I. -1>. 

-J' "J' J« J* J« r 1 - r 1 ' }>• <-i- - 1 ' i. l i. -1. i. -1). {-i. -i. 1. 1. 1. -1. 1. 1}. 

-1, -1, 1, 1, 1, 1, -1, -1}, {-1, -1, x, 1, 1, 1, -1, 1}, {-1, -x, 1, l. 1, 1, 1, -X}, 

"I' T 1 ' !• 1 ' 1 ' !' LI}. JC,- 1 ', 1 ' _1 * _1 * X ' - 1 ' - 1 ' t" 1 ' - 1 ' "I- 1' "I- -1. 1). 

-1, 1, -1, -1, 1, -1, 1, -1}, {-1, x, -1, -1, 1, -1, 1, 1}, 

-1, 1, -1, -1, 1, 1, -1, -1}, {-1, 1, _x, -X, 1, 1, -1, 1}, {-1, 1, -1, -x, 1, 1, 1, -1}, 

-\- }• -}• i. i. J- 1 >; {-i. i. -i. i. i. -i. -i. -i). {-i. i. -i. i. i, -1. -1, 1}, 

-1, 1, -1, 1, 1, -1, 1, -1}, {-1, X, -1, 1, 1, -1, 1, 1}, {-X, 1, -1, 1, 1, 1, -X, -1}, 

V 7 1 ' J* }• x ' - 1 ; i). ,{-i. i. -i. i. i. i. i. -i). (-1. i. -i. i. l. l. i, i}, 

-1, 1. 1. -1. 1, -l. -l, -1}, {-1, 1, 1, -1, l. -l, _x, 1}, {-1, 1, 1, -x, 1, -i/i, -1), 
-1. 1, 1, -1, 1, -1, 1, 1}, {-X, 1, 1, -1, 1, 1, -1, -1}, X, 1, -1, 1/ X, -1, X}, 

}' }' r 1 ', 1 - }' 1 ' i» L -i. i. i. i. i). {-i. i. i. i. i. -i. -l. -i>. 

-1, 1, 1, 1, 1. -1. -1, 1}, {-1, X, 1, 1, 1, -1, 1, -1}, {-1, x, 1, 1, 1, -1, 1, 1}, 

-1, 1. 1, 1, 1, 1, -1, -1}, {-1, 1, 1, 1, i, 1, -x, 1}, {-X, 1, 1, 1, 1, 1, 1, -1}. 

-1, 1. 1, 1, 1, 1, 1, 1}, {1, -1, -1, -1, 1, -X, -1, _x } , {1 , -X, -1, -1, 1, -1, -1, 1}, 

1, -1, -1, -1, 1, -1, 1, -1}, {1, -1, -1, -1, X, -1, 1, 1}, " 

-J* -J- -J- J- i' r 1 ' {1 - _1 ' _1 ' - 1 ' x ' x ' 1} - {1 - -i- -i- -I- i< !' i' -i}. 
J- "J* r 1 ', 1 - J« J* 1 >', <i- -i' -L i. i.'-i. -i. -1). {i. -i. -i. i. i. -i. -i. i>. 

1, -1, -1, 1, 1, -1, 1, -1}, {l. -l, -1, 1, 1, -1, x. 1}, {1, -X, -1, 1, 1, 1, -1, -1}, 

1, -1, -1, 1, 1, 1, -1, 1}, {1, -1, -1, 1, x, 1, 1, -1). {1, -1, -x, 1, 1, 1, 1, X}, 

1. -1, 1, -1, 1, -1. -1, -1}, {1, -1, 1, -1, 1, -1, -1, 1), {1 , -1, 1, -X, 1, -1,1. -1}, 

1, -1, 1, -1, 1, -1, 1, 1), {1, -1, 1, -1, 1, x, -1, -1}, {1, -1, 1, -1, 1, x, -1, 1), ' 

1, -1, 1, -1, 1, 1, 1, -1}, {1, -1, 1, -1, x, 1, 1, 1}, {1, -1, 1, 1, x, -1, -1. -1}, 

1, -1. 1, 1. 1. -1, -1, 1}. {1, -1, 1, 1, 1. -1, X, -1}, {X, -1, 1, 1, 1, -1, 1, 1), 

1, -1. 1, 1, 1. 1, -1. -1}, {1, -1, 1, 1, i; 1, -X, 1}. {X, -1. 1, 1, 1, 1, X, -1>, 

1, -1, 1, 1, 1, 1, 1, 1}, {1, 1, -1, -1, 1, -1, _x, -X}, {1, 1, -1, -1, 1, -1, -x, 1}, 

1, 1, -1, -1, 1, -1, 1, -1}, {1, 1. -1, -1, 1, -x, 1, 1}, {1, 1, -X, -1, 1, 1, -X, -1}, 

1, 1, -1, -1, 1, 1, -1, 1}, {1, 1, -1, -1, 1, x, 1, -1}, {1, 1, -1, -x, 1, 1, 1, X}, 

1, 1, -1, 1. 1, -1, -1, -1}, {1, 1, -1, 1, 1, -x, -1, 1}, {1, 1, -x, 1, 1, -1, 1, il), 

1, 1, -1, 1, 1, -1, 1. 1}, {1, 1, -X, 1, 1, 1, -1. -1}, {1, 1, -1, x, 1, 1, -1, 1}, 

1, 1, -1, 1, 1, 1, 1, -1}, {1, 1, -1, 1, 1, 1, 1, X}, {1, 1, 1, -1, 1, -1, -x, -X}, 

}- }■ i- -i' i- r 1 ' -i- 1 >- fi' "i. 1 - -i' 1. -i. i. -1). {i. i. i. -i. i. -i. 1. i}. 

1, 1, X, -1, 1, 1, -1, -1}. {1, 1. 1, -1, 1, 1, -x, 1}, {X, 1, 1, -1, 1, 1, X, -1}, 

1. 1. 1, -1, 1, 1, 1, 1}, {1, 1, 1, 1, 1, -X, -1, -1}, {1, X, 1. 1, 1, -1, -1, X}, 

1, 1, 1, 1, 1, -1, 1. -1}, {1, 1, 1, 1, 1, -X, 1, 1}, {1, 1, X, 1, 1, 1, -1, -1}, 

1, 1, 1, 1, 1, 1, -1, 1}, {1, 1, 1, 1, X, 1, 1, -1}, {1, l, 1, 1, X, 1, 1, 1}} 
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xl = Map[ModulationValue[FiltPulse[8] ] [#] [0.5 T]&, 

^Select [Table [i, {i, 1, 2*8}] // Map [ConvertToBitSeq[8] , #]&, (# [[5]] == 1)&] ] 



{0.749266 
0.829796 
.776885 
.802178 
.759521 
.819541 
.78714, 
0.791922 
0.749266 
0.829796 
0.776885 
0.802178, 
0.759521 
0.819541, 
0.78714, 
0.791922, 



74922, 
82975, 
776839 
802131 
759475 
819495 
0.787094, 
0.791876 
74922, 
82975, 
776839 
802131 
759475 
0.819495 
0.787094, 
0.791876 



Map [1/2 - a/2 Erf 



0.749173, 
0.829703, 
, 0.776791 
, 0.802084 
, 0.759428 
, 0.819447 
0.787047, 
, 0/791829 
0.749173, 
0.829703, 
, 0.776791 
, 0.802084 
, 0.759428 
, 0.819447 
0.787047, 
, 0.791829 



0.749219, 
0 . 829749, 
, 0.776838 
, 0.80213, 
, 0.759474 
, 0.819494 
0.787093, 
, 0.791875 
0.749219, 
0.829749, 
. 0.776838 
, 0.80213, 
, 0.759474 
, 0.819494 
0.787093, 
, 0.791875 



0.711482, 

0.792012, 
, 0.739101 

0.764394, 
, 0.721738 
, 0.781757 

0.749356, 
, 0.754138 

0.711482, 

0.792012, 
, 0.739101 

0.764394, 
, 0.721738 
, 0.781757 

0.749356, 
, 0.754138 



0.711529, 0 

0.792059, 0 
, 0.739147, 

0.76444, 0. 
, 0.721784, 
, 0.781804, 

0.749403, 0 
, 0.754185, 

0.711529, 0 

0.792059, 0 
, 0.739147, 

0.76444, 0 
, 0.721784, 
, 0.781804, 

0.749403, 0 
, 0.754185, 



711482, 0 
792012, 0 
0.7391, 0. 
764393, 0. 
0.721737, 
0.781756, 

749355, 0 
0.754138, 
.711482, 0 
.792012, 0 
0.7391, 0. 
764393, 0. 
0.721737, 
0.781756, 
.749355, 0 
0.754138, 



711435, 
791965, 
739054, 
764347, 
0.72169, 
0.78171, 
749309, 
0.754091, 
.711435, 
.791965, 
739054, 
764347, 
0.72169, 
0.78171, 
749309, 
0.754091} 



// Apply [Plus, #]& // #/128 &; 



Ber[a_] : = Evaluate [%4 8 ] 
Ber [0.01] 

0.495 

Erf [Infinity] 

1 



D[Erf [x], x] 

2 E' x2 



? Erf 



Erf[z] gives the error function erf(z). Erf[z0, zl] gives the generalized error 
function erf (zl) - erf(zO). 

Plot [Erf [x] , {x, 0, 10}] 



0.9999 
0.9998 
0.9997 
0.9996 
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Graphics - 



TFiltPulse 



Global 'FiltPulse 
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T := 3 / 812500 

Null 

?? T 

This is the symbol period 
?? T 

This is the symbol period 
T := 3 / 812500 

Clear[T, xO, xl, x2, x3, x4, x5, x6, x7] 

j(a H . 4 *a KO *a B . 2 *a„.j ♦ a N ) p^^o] [ £ T j + jfa„^ + a H . a ♦ « s _ 2 ) ^^^j [ 2 T + <5T] + 

J**-* Pulse[0] [4T + <5T] + j<-**-s) p u lse[0][6T + <5T] + 

j(-a M . 5 -a«. 8 -a s . 7 ) p^isefoHST + ST] + tf<«»-4 + a«.» ♦ puise[l][T + ST] + 

j(«no> Pulse[l] [3T + ST] + J<-*»-«> Pulse[l] [5 T + <5T] / . { a„ -* xO, - xl, 

a N . 2 -+ x2, a N - 3 -*> x3, a N _ 4 x4, a H _ 5 x5, a N _ 6 -+ x6, a N _ 7 x7} 

(1. I)* 0 ** 1 -* 2 - 3 *** Pulse[0] [<5T] + (1. i)^x 3+ x4 Pulse[0] [2T + <5T] + 
(1. I) x4 Pulse[0] [4T + <5T] + (1. I)" x5 Pulse[0] [6T + «5T] + (1. !)-x5-x6-x7 Pu i se[C j] [8T + {5T] 
(1. i)*i+x3+x4 Pulse [1] [T + <5T] + (1. I) x3 Pulse [1] [3T + 5T] + (1. I)-* 6 Pulse [1] [5 T + «5T] 

*<««-»> Pulse [1] [3 T + <5T] /. {a N -+ xO, a,,.! -* xl, a N _ 2 x2, a N _ 3 -+ x3, 
a H . 4 -► x4, a N _ 5 -+ x5, a N _ 6 -+ x6, a N _ 7 -+ x7} 

(1. I) x3 Pulse [1] [3T + c5T] 

(1. * I)* 3 Pulsefl] [3 T + ST] 

(1. I) x3 Pulse [1] [3 T + <5T] 
? J 

J = € i8 
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BeginPackage [ "LaurentFunc tions x " ] 

T: : usage = "This is the symbol period" 

BT::usage = "This is the usual product" 

h: :usage = "This is the raw gaussian pulse" 

$::usaga = "This denotes modulation index Pi" 

ip: : us age ="^[L,t] is Laurents function" 

hFiltered: rusags = "This is the filtered gaussian pulse" 

PhaseAngle: : usage = "This function takes some time to calculate. The following 

code will draw a graph of the functionPhaseFunction[t__] =N[ £ L t ] ; 
phasepoints =Table [ {t , PhaseFunction [t T] } //N, {t, 0, L ,1/40}]; 
ListPlot [phasepoints. Plot Joined -> True]" 



S : 


: : usage = 


»s 


= Sin[B] " 


J: 


: us age = 


"J 




C: 


: usage = 


H C 


= Cos [$] " 


M: 


rusage = 


"M 


= 2 L " 1 » 



ModulationXndex: : usage = "Modulationlndex = h" 

LaurentS: :usage = «LaurentS[L] [n] [t] = Sin[ tfr[L,t + n T] ] / S" 

LaurentC: :usaga = "LaurantC [L] [K] [t] is Laurent s C K , t " 

AlphaKI: rusage = "AlphaKI [LL] [K, i] is Laurents a Kt i» 

LaurentLK: rusage = "LaurentLK[L] [K] gives the support of C K , t " 

The start of Modulator Definitions 

ANKInitialStateSetUp: :usage = 
"ANKInitialStateSetUp[I,] [K] [InitBitSeq, AccumulatedPhase] sets up the sequence 
of prior states of A K , H . that the modulator went thgough to get to the 
constelation point specified by AccumulatedPhase which is really A 0 , 0 " 

AKN: rusage = "AKN[L] [K] [ {State, AccumulatedPhase} ] defines A K , H in terms of A 0 , K " 

ModulatingPulse: : usage = 
"The Pulse is assumed to have the following structure Pulse [L] [K] [t] « 

NumberOf Curves: -.usage = "The number of pulses used by the modulator" 

Samplinglnterval: : usage = 
"Samplinglnterval is the interval between samples of the output of the modulator" 

InitialState: :usage = "InitialState is the set of bits that are assumed 
to be present before i.e {a. 1/ a_ 2 , •••}" 

StartingQuadrant: rusage = "StartingQuadrant = A 0 ,_ ia nd is a number" 

Modulator: rusage = "Modulator [L] [BitSeq, Opts] assumes the following 
default options StartingQuadrant -> 0, InitialState -> Table[l,{i, 
1/ 20} ], Samplinglnterval T/3 2 , NumberOf Curves -> 4 , ModulatingPulse-* 
LaurentC . The Pulse is assumed to have the following structure Pulse [L] [K] [t] « 

Start of the Receiver Functions 
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FiltPulse: :usage = "The default pulse and is called by FiltPulse[L] [K] [t] ■ 

Sync Sample :: us age = "Given that the sampling interval is T/32,then sync 

sample has rang -16 to 16 and this moves the point used to demodulate " 

Receiver: :usage = '» Receiver [L] [Input Seg f Opts] assumes the following- 
default options {StartingQuadrant-*0 # InitialState-» { 1, 1, 1, 1, 1, 1, 1, 1, 
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1} , Samplinglnterval -> T/32 , KodulatingPulse -+ 
FiltPulse, SyncSrjsiple -> 0 , K\LnberOf Curves -> 2 " 

ReceiverProper: :usage = "Receiver Prope r [ L] [ S tart ingQuadr ant, InitialState, 
Sampling Interval , Modul a t ingPul s e , Sync S ample , NumberO f Curve s , Input Seq] 
is called by Receiver after all tha options have been resolved" 

Eogin [ " * Private * " ] 

V^og[2] 

a : = 

2 7TBT 



h[t_] := 



'2 7T a T 



Ideally we would and did define the effect of the convolution of h[t] by the formula below. However* this version of 
Mathematica gives. an error. 

f h[t- r] 



kf±i [t__] := Release [Module [ {r} , J ' ^ — ~ dt]] 



hFiltered[PulseWidthJ [t_] := Module [ {xl # x2, x3}, xl = 

. r rr fT h[tl-t] ,- r PulseWidth PulseWidth T 
40] & )[Table[{tl, J _L_J- d r}, {tl, , JL}]], 

x2 = Interpolation [xl] ; x2 [t] ] 

B : = N [Modv^lat ionlndex tt] 

C := Ccs[5] ; 
S : = Sin [£] ; 

J := ( <e j * // Chop) ; 
M := 2 1 *- 1 ; 

PhaseAngle[L__] [t_] /; t S 0 : = 0 
PhaseAngle[L_J [t_] /;tiLT:= 3 

PhaseAngle [L_] [t_] := 

PhaseAngle[L] [t_] = Module [{xl, x2, x3,,x4, x5, x6}, xl = hFiltered [3 I* T] [ tl - -^-1; 

2 J 

x2 = Table [{t2, $ Evaluate [xl] dltl}, {t2, 0, LT, }]; Interpolation [x2] [t]l 

J-LT 100 

^[L_, t_] /; 0< t < LT := Phase Angle [L] [t] 

^[L_, t_] /; 2LT > t * LT : = 9 - PhaseAngle [L] [ t - L T] 

We need to put this to avoid the function being extrapolated 

</r[L_ , t_] /; ! (0 < t < LT) && ! (2 LT > t * LT) := 0 
LaurentS[L_] [n_] [t_] := Sin[^[L, t + nT]]/S 
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AlphaKI [LL_] [K_, i_] /; (0 < i < LL) && (0 <= K < 2"" 1 ) : = 

r ( .KNum - x2 \ 
Module [ x2, x3, KNum} , xl : = ^x2 = Mod [KNum, 2] ; KNum = ; x2 I ; 

KNum s K; x3 = Table [xl, {ii, 0 , LL - 1} ] ; x3 [[i]I ] 

LaurentLK [L_] [K_] ;= 
Module[{xl}, xl = Table [L (2 - AlphaKI [L] [K, ii] ) - ii, {ii, 1, L - 1} ] ; Min [xl] ] 

LaurentC[L_] [K_] [t_] /; 0 <= K < 2 L : = 

L-l 

UurentS[L] [0] [t] J~J LaurentG[L] [ii + L AlphaKI [L] [K, ii] ] [t] 

£i=l 

The start of the modulator function 

ANKInitialStateSetUp[L_] [K_] [InitBitSeq_, AccumulatedPhaseJ : = 
Module [{xl, x2, x3, x4, x5, acuphase, initbitseq}, 
initbitseq = InitBitSeg; 
acuphase = AccumulatedPhase; ■ 
UpdateSeg : = 

Module [{}, xl = acuphase - Sum [initbitseq [ [i ] ] AlphaKI [L] [K, i] , {i, l, l -1}]; 
acuphase = acuphase - First [initbitseq] ; initbitseq = Rest[ initbitseq]; xl] ; 
Table [ UpdateSeq, {i, 1, Laurent LK [L] [K] } ] ] 

AKN[L_] [K_] [{State_, Aoc-;mulatedPhase_} ] : = 
AccumulatedPhase - Sum [State [ [i + 1] ] AlphaKZ [L] [K, i] , {i, 1, L -1}] 

Options [Modulator] : = {StartingQuadrant -> 0, InitialState Table [1, {i, 1, 20}], 
Saiaplinglnterval -> T/32, ItoiberOfCurves -+ 4, ModulatingPulse-* LaurentC) 

Modulator [L_] [BitSeq_, Opts ] : = 

Module [ 

{xl, x2, x3, x4, x5, x6, state, AccumalatedPhase, seq, AKKState, Curves, Pulse}, 
xl = Sainplinglnterval / . {Opts} / . Options [Modulator] ; 
state = InitialState /. {Opts} /. Options [Modulator] ; 
x3 = StartingQuadrant/. {Opts} /. Options [Modulator] ; 
x4 = Sair.pl inglnterval /•. {Opts} /. Options [Modulator] ; 
Pulse = ModulatingPulse/. {Opts} /. Options [Modulator] ; 
Curves = (Ku^erOf Curves / . {Opts} / - Options [Modulator] ) -1; 
AccumulatedPhase = x3 ; 
seq ~ BitSeq; 
Table [ 

AKNState[K] = ANKInitialStateSetUp [L] [K] [state, AccumulatedPhase], {K, 0, Curves}]; 
x5 : = Module[{}, state = Join[ {First [seq] } , Drop[state, -1]]; 
AccumulatedPhase = AccumulatedPhase + First [seq]/ 
seq = Rest [seq]; 

Table [AKNState [K] = Join [ {AKN [L] [K] [{state, AccumulatedPhase}]}, 

Drop[AKNState[K] , -1]], {K, 0, Curves}); 
x6[r_] = Sum[Sum[ ( j)*™tat. [k] [ [i ♦ i] j PulsetL] [K] [r + ± T] f 

{i, 0, LaurentLK[L] [K] -1}], {K, 0, Curves}]; 
Table [x6[r], {r, 0, T - x4, x4}]]; 
Table [x5, {kk , 1, Length [BitSeq] } ] //Flatten] 

Options [Receiver] := {StartingQuadrant -» 0, 

InitialState-* {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1}, Saroplinglnteir< 
ModulatingPulse -> FiltPulse, SyncSample -> 0, NumberOf Curves -> 2}; 

Receiver [L__] [Input Seq_, Opts ] := v 

Module [{xl, x2, x3, x4, x5, x6} , 

xl = StartingQuadrant/. {Opts} / . Options [Receiver] ; 
x2 = InitialState/. {Opts} /. Options [Receiver] ; 
x3 = Samplinglnterval / . {Opts} /. Options [Receiver] ; 
x4 = ModulatingPulse/. {Opts} /. Options [Receiver] ; 
x5 = SyncSample/. {Opts} /. Options [Receiver] ; 
x6 = NumberOf Curves / . {Opts} /. Options [Receiver] ; 
ReceiverProper [L] [xl, x2, x3, x4, x5, x6, Input Seq] ] 
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ReceiverProper[L_] [ StartingQuadrant„, InitialState_, Samplinglnterval_, 
Modulating?ulse_, SyncSample_, NimberOfCurves_, InputSeq_] : = 

Module [{xl, sgn, ReceivefiSeq, ExpectedValue, seg, ReceiveNext, ID, Si), 
xl = T / Sainplinglnterval; 

ReceivedSeg = Partition [Input Seq, xl] // Transpose //#[ [Sync Sample + 1]]& 
ExpectedValua = 

Modulatir.g?ulse[L) [0] [ (LaurentLK[L] [0] / 2) T + Sync Saiaple Sampling Inters 
& = StartincQuadrant; 
sgn = 0; 
ID = {}; 

seq = ReceivedSeq; ^ 
ReceiveNext : = 
Module [{xl, x2>, 

xl = ( (- 1) SffU J 91 First [seq]) // Ira; 
If[AJDs[xl - ExpectedValue ] <= Xbs [xl + ExpectedValue ] , 

3D = Join[ Z>/ {!>]; # = + 1, ID = .Toin[ 3D, {-1}]; & = & - 1 ] ; 
seq = Rest [seq] ; sgn = Modfsgn + 1, 2]]; 
Table [ReceiveNext, {i, 1, Length [ReceivedSeq] }] ; 
3D] 

End [ ] 



EndPackage [ ] 



Lauren tReceive rWithLook UP.). 



1 



Needs [ "LaurentFunctions * " ] 



RuleDelayed: :rhs : Pattern t_ appears on the right-hand side of rule 

PhaseAngle[L_] [t_] (PhaseAngle [L] [t_] = Module [{xl, x2, x3, x4 , x5, x6}, «1»]) . 

Needs [ "LaurentNotationTest 1 " ] 

Needs: :nocont : Context LaurentNotationTest* was not created when Needs was evaluated. 



Information on the functions used can t?_e. obtained using help. 



Names [ N LaurentFunctions ^ * N ] 

{AKN, AlphaKI, ANKInitialStateSetUp, BT, FiltPulse, h, hFiltered, InitialState, J, 
LaurentC, LaurentLK, LaurentS, M, ModulatingPulse, Modulationlndex, Modulator, 
Number Of Curves, PhaseAngle, Receiver, ReceiverProper, S, Samp ling Interval , 
StartingQuadrant, SyncSample, T, C, 

3 

T := 

812500 
BT := 0.3 



1 

Modulationlndex : = — 
2 



<< ModulatorData .m; 



RandomBitSeq 

{1, 1, -1, -1, -1, 1, 1, -1, 1, -1, 1, .1, 1, -1, -1, -i, _i, i, i, i, _ 1/ _ X/ lf ! x 

-1, 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, -1, 1, -1, -1, i, _i, lt lf lt _ ± 

1, -1, 1, -1, 1, -1, 1, 1, -1, -1, 1, -1, 1, -1, 1, 1, i, i, -i, lf 
1, 1, 1, -1, 1, -1, -1, 1, 1, 1, -1, -1, 1, l, l, l, l, i, i, i, i, _i f lt -1, 1, -1} 
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With L = 8, and M = 2 we can utilise only the main values. It so happens that when K = 1, the dominant value occurs at 
2.5T in the filtered pulse. The three dominant values when K =0 occurs at 3.5T, 4.5T and 5.5T. The value of the pulse 
at 6.5T is smaller than the value of the pulse with K = 1 at 2.5 T. Thus the expression becomes 



S N T + AT = > J] jAK " N " n " C *.n'T + AT 

§nt + at = J* 0 -"-* Pulse[0][3T + ST] + Pulse[0] [4 T + <5T] + 

J X °' N - 5 Pulse [0] [5T + 6T] + j*i.B-a p u l S e [1] [2 T + <5T] 

Siit^at = J* 0 - 11 - 5 Pulse[0)[3T + 6TJ + Pu i se [0] [4 T + <5T] + 

j* 0 ,m-5 Pulse[0]"[5T + 5T] + j^-M-s^K.,) Pulse[l][2T + <5T] 

Given that a N _ 3 has been already decided, we can precalculate the possibilities and store them. 

S N t*at = J* 0 '"- 5 ( j{°«-** a «-3> PulsetO][3T + <5T] + CT 0 *-' Pulse[0J[4T + <5T] + 
Pulse [0] [5 T + <5T] + Ji°s- «••»-» • °n-2*°n- s > Pulse[l][2T + <5T] ) 

But J 2 =1 and so we get 

S HStM B ( jtas-i^aM.j) Pu i se [o][3T + «5T] + CT 3 "-* Pulse[0][4T + <5T] + 

Pulse [0] [5 T + <5T] + ji°h-** ) Pulss[l][2T + <5T] ) 

LookUpTable [Pulse_, <5T_] : = 
Module [{xl, x2, x3 # x4), 
xl = {{-1, -1, -1}, {-1, -1, 1}, 

{-1, 1, -1), {-1, 1, 1}, {1, -1, -1}, {1, -1, 1}, {1, 1, -l}, {i, i, 1}}; 
x2 [bitsecr ] := 

( J bit«.q[C3]]*b±t«,[[an p u lse[03[3T + <5T] + jbito^t m ] Pulse[0][4T + + 
Pulse[0][5T + <5T] + j* 1 * 8 *^ [*] ] ♦ bitseq[ W 3 Pulse[l][2T + <5T] ) ; 
{xl, Map[x2[#]&, xl]} //Transpose] 

tab = LookUpTable [FiltPulse [-8] , 0 . 5 T] 

{{{-1, -1, -1}, -0. 00122183 - 0.770616 1}, {{-1, -1, 1}, 0.702339 + 0.770616 1}, 
{{-1, 1, -1}, 0.614125 - 0.770616 1}, {{-1, 1, 1}, 0.0869922 + 0.770616 1}, 
{{1, -1, -1}, 0.0869922 - 0.770616 1}, {{1, -1, 1}, 0.614125 + 0.770616 1}, 
{{1, 1, -1}, 0.702339 - 0.770616 1}, {{1, 1, 1}, -0.00122183 + 0.770616 1}} 

Sort[tab, 3]] >£2[[1, 3]]&] 

{.{{1, 1, 1}, -0.00122183 + 0.770616 1}, {{1, -1, 1}, 0.614125 + 0.770616 1}, 
{{-1, 1, 1}, 0.0869922 + 0.770616 1}, {{-1, -1, 1}, 0.7 02339 + 0.770616 1}, 
{{1, 1, -1}, 0.702339 - 0.770616 1}, {{1, -1, -1}, 0.0869922 - 0.770616 1}, 
{{-1, 1, -1}, 0.614125 - 0.770616 1}, {{-1, -1, -1}, -0.00122183 - 0.770616 1}} 

FiltPulse [8] [0] [4 T + 0 . 5 T] 

0.770616 



Now we build a receiver! ! First generate the modulated sequence 
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Options [Modulator]. 

[startingQuadrant -> 0, 
InitialState-* {1, 1, 1, 1, 1, 1, 1, 1, l, l, l, l, i ( i, 1# 1# lf 1} ( 

Samplinglnterval -» 260 q Q00Q , NumberOf Curves -* 4 , ModulatingPulse -* Laurentc} 



26000000 



1 

32 



tom2 = Modulator [Ii] [Random3itSeg, NumberOf Curves 2, ModulatingPulse FiltPulse] ; 
ListPlot[{Re[tom2], Im[tom2]) //Transpose, PlotJoined -> True] 

- Graphics - 
RandomBitSeq 

{1, 1, -1, -1, -1, 1, 1, -1, 1, -1, 1, -i, i, _l, _i 7 -l, -i, i, !, lt _ 1/ _ a _ 1 lf ! ! 

r 1 - r 1 * :V - 1 * "V V' J- - 1 * ^ - 1 * -i. -i, i. 1/1. i. -1. 

1, 1, 1, -1, 1, -1, -1, 1, 1, i f -1, -1, i, i, i, i, !, lf 1§ 1# lt „ 1# _ lt Xt _ lt lf _ 1} 
Receiver [L] [tom2, StartingQuadrant -> 4] 

{1, 1, 1, 1, 1, 1, -1, -1, -1, 1, l, -1, 1, -l, l, _i, i, „i, . x 1# x 3L ^ _ 1 

-1*1.1, 1. -1, 1. -1. "I, 1, -1, -1,1, 1, 1, 1, -1,1, -1, 1, -1, 1,-1 -i, i, ' ' 

x \ 1 s 1 'r lf , 1 'n" 1 i 2 ; x ' x ' 2 ' - 1 ' - 1 ' x - ^ x < ^ i- 1* -i' 

-i, -i, i, -1, i, i f i, -i, i, -i, -i, i, i, i, -i, -i, x, i, !, lf lt lt lt lf lf _ 1# _ 1} 

We now calculate the BER given the and L = 8. We write all the possible terms 

• j(«».4*««-»*««- a * a M .i) pul se [0] [T + <5T] + jf fl »-**«H-3 * pui se [0] [2 T + <5T] + 

j(« w -**ok-3) pulse [0] [3 T + 6T] + J**-* Pulse[0] [4 T + <5T] + Pulse[0] [5T + <5T] + 
j(-a»-s) Pulse[0][6T + <5T] + j<-°n- 5 - a„. e ) Pulse [0][7T + <5T] + 

j(-a N . s -a N . fi -a H _ 7 ) p u l se [0] [8 T + <5T] + ♦ °n - s ♦ «n - i ♦ »r> - °n - i > pul se [1J [ <5t] + 

tf<«H-4*««0 ♦ ♦««-! - ) Pulse[l] [T + <5T] + ,J(«N-4*ON.3 * °W - 2 * a M - J ) p U lse[l] [2 T + ST] + 

j£°n-**o»-3 > puise[l] [3 T + <5T] + j(°n-4- o h-s ) Pulse [1] [4 T + <5T] + 
j(-°h-« ) Pulse [1] [5 T + <5T] + j(-«h- s * °h - 7 ) Pulse [1] [6 T + <5T] ) 

We select the imaginary 

jAo. H _s ( J^-^ 0 "-^ 0 "^ +a »-i Pulse[0] [ <5T] + j{°».4*c N . 3+ a K . 3 ) Pulse[0]E2T + ^ + 

' * Pulse[0][4T + <5T] + j(-«n-s) p u lse[0][6T + <5T] + 

Cr(««-A*OH-3 ) pulse [1] [3 T + <5T] + Jf-'N-* ) Pulse[l][5T + <5T] ) 

ModulationValue[Pulse_] [{x0_, xl_, ac2_, x3_, x4_, y.5_, x6., x7_) ] [<5T_] i = 
(1.* xjxo.xi~2.x3~4 pulse [0] [5T] + I) x^x3.x* Pulse[0] [2 T + 5T] + 

(1. * I) x * Pulse[0] [4 T + <5T] + (1. * I) ~ x5 Pulse[0] [€ T + 5T] + 

JJ-S-xS-xT p ulse[0] r 8 T + 6T] + (1^ T)^~3.x 4 p ulse[1] [T+<5T] + 

I) 543 Pulse [1] [3T + 6T] + (1.- I) ~ xS Pulse [1] [5T+<5T] 
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MoaulationValue[FiltPulse[8J] £{-1, -1, -1, -1, _i, _ X/ _ 1# -1)][0.5T] // N 
-0.754091 1 



ConvertToBitSeq[Depth_] [Num_] := 
Module [{bit, n} , 
n = Num; 

ConvertNext := (bit = Mod[n, 2]; n = Floor[n/2]; bit); 
( Table [ConvertNext, {i, 1, Depth}] //Reverse) /. {0 -> -1)]; 

ConvertToBitSeg[8] [2] *" * "" 

{-1, -1, -1, -1, -1, -1, 1, -1} 



Table[i, {i, 1, 2 A 8 } ] // Map [ConvertToBitSeq [8 ] , #]&; 

Select [Table [i, {i, 1, 2 A 8}] // Map [ConvertToBitSeq [8] , #]&, (#[[5]] 1) &] 



-1, -1, -1, -1, 1, -1, -1, -1}, 

-1, -i, -1, -i, i, -l, -i, i), {-i, -i, -i, -i, i, -1,-1, -i}, 

-1, -1, -1, -1, 1, -1, 1, 1}, {-i, -i, -l, -i, i, !, 
-1, -1, -1, -1, 1, 1, -1, 1}, {-l, -i, -i, -i, 1# X/ 
-1, -1, -1, -1, 1/1, 1, 1), {-l, -l, -i, i, i, -i, 
-1, -1, -1, 1, 1, -1, -1, 1), {-l, -i, -i # i, a, !, 

-1, -1, -1, 1, 1, -1, 1, 1}, {-1, -l, -i, i; i, _ lf { _ lt _ lt _ lt 1 J, x _ x 1} 

-i, -l, i, -l, i # -i, -i, l), {-i, -i, i, -i, i, -i}, " 

-1, -1, 1, -1, 1, -1, 1, 1}, {-1, -1, 1, -1, 1, 1, -i, -i}, _i ( i t _ 1# lf x .J 1} 

T 1 ' lt 1 ' l ' t" 1 - X > X > X ' X >' f" 1 ' a ' ^ -i. -i}. 

-l, -i, i, i, i, -l, -l, l}, {-i, i, i, -i, i, -l), {-i, -i, i, i, _ 1; lt 1} ; 

-1, -1, 1, 1, 1, 1, -1, -1), {-1, -1, 1, l, l, _i, i}/ , lt lt lt lt lf ! _ 1} 

-1, -1, 1, 1, 1, 1, 1, 1}, {-1, 1, -1, -1, 1, -1, -l, -ij, {-!, i, . lf lt .i # _{ 1} 

-1, 1, -1, -1, 1, -1, 1, -1), {-l, l, -i, -i, i, i, 1} , 

-1, 1, -1, -1, 1, 1, -1, -1}, {-1, 1, -1, -1, i, i, -i, i }/ lt _! _ 1/ x x x „ 1} 

-i, i, -i, i, i,. -i f i, -i}, {-i, i, -i, i, i, -i, i, i} # {_i, i, -i, i, i, i, 

-1, 1, -1, 1, 1, 1, -1, 1}, {-1, 1, -1, 1, 1; 1, 1, -1), {-1, 1, -1, 1, 1, lf 1# 1} 

-1, 1, 1, -1, 1, -1, -1, -1}, {-l, 1, 1, -1, 1, -l, -i, i } , {_!, 1( lf . 1# lf _ 1 / lf „ 1} 

-1, 1, 1, -1, 1, -1, 1, 1}, {-1, 1, 1, -i, 1, 1,-1, -l}, {-l, i, -i, i, i, .! 1} 

-1, 1, 1, -1, 1, 1. 1, -1), {-1, 1, 1, -1, i, i, i, i}/ X/ 1# lf lt -]L . 1 J 1} ' 

-1, 1, 1, 1, 1, -1, -1, 1}, {-1, 1, 1, 1, l, -l, _!}, {-!, i/ 1# lf lt 1} 

-1, 1, 1, 1, 1, 1, -1, -1}, {-1, 1, l # l, -i, 1} , ! 1# x x x ! „ 1} 

-l, i, i # i, i, i, i, i}, {i, -i, -i, -i, i, -i, -i, { i, _i, _ 1( . lr 1§ 1}j 

1, -1, -1, -1, 1, -1, 1, -1}, {l, -1, -l, -l # l, -i, i}, 

1, -1, -1, -1, 1, 1, -1, -1}, {1, -1, -1, -l, l, l, -i, i}, { i, _! lf x x _ 1} 

i, -l, i, i, i, i}, {l. -l, -l, i, i, -i, -i, -i), {i/ -i, -i, i, x, 1} ; 

l f -1, -1, 1, 1, -1, 1, -1), {1, -1, -l, l, i, -i, i, i } , {lg _ lf _ lt lt l, 1, -i, 

1, -1, -1, 1, 1, 1, -1, 1}, {l, -1, -l # i, i # 1, 1, {lf _i, lf 1# 1 x 1} 

1, -1, 1, -1, 1, -1, -1, -1}, {l # -1, 1, -l, l, -i, -i, i } , { i, -i, !, _ lt lf _! lt . 1} 

1, -1, 1, -1, 1, -1, 1, 1}, {1, -1, 1, -1, 1, 1, -1, {li _i, X , _ X/ x ! _ lf 1} 

i, -l, i, -l, i, i, i, -l), {i, -l, i, -l, i, i, i}, { i, i, x, 1/ . 1# 

1, -1, 1, 1, 1, -1, -1, 1}, {1, -1, 1, 1, l, -1, 1, (i, -i # i, i, lt _ 1# lt 1} 

i, -l, i, i, i, i, -i, -i}, {i, -l, i, i, i, i, -l, i}, { i, -i, i, i, !, x, lf 

1, -1, 1, 1, 1, 1, 1, 1}, {1, 1, -1, -l, -i, -i, { i, _ 1# _ lt lt 1} 

1, 1, -1, -1, 1, -1, 1, -1}, {1, 1, -l, -l, 1, -i, i, 1} , {1 , lt _ 1/ _ lt lt lt _ lf _ 1} 

1, 1, -1, -1, 1, 1, -1, 1}, {1, 1, -1, -1, i, i, i, { i, lf . lf _ lt lt x x 1} 

1, 1, -1, 1, 1, -1, -1, -1}, {1, 1, -l, l, i, -i, 1} , { i, Xt _ Xf lt lt „ X/ ! i 1} 

1, 1, -1, 1, 1, -1, 1, 1 }/ {1, 1, -l, l, 1, 1, -i, {1/ lt _ lf lf 1# lf - lf l ' )f 

1, 1, -1, 1, 1, 1, 1, -1} # {1, 1, -1, 1, 1, 1, i), {1 , X/ i, _i ( 
1, 1, 1, -1, l f -1, -1, 1}, {1, 1, l, -l, i, -i, i, {1 , lt lf .! 2 _ lf lf 1} 

i, i f i, -l, i, i, -l, -l}, {i, i, i, -i, i, i, _i, 1} , {1 , lt lt _ 2/ 1# 1# lt 

1, 1, 1, -1, 1, 1, 1, 1}, {1, l, 1, 1, 1, _i # -i, {1/ i, lt lt lt _ 1( _ lf 1} 

1, 1, 1, 1, 1, -1, 1, -1), {1, 1, 1, 1, -l, i, i } , { i, i, i, i, !, lf _ 1} 
1, 1, 1, 1, 1, 1, -1, 1}, {1, 1, 1, 1, 1, 1, 1, -1}, {i, i, i, i, i, i, 1/ 1)} 
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xl = Map[ModulationValue[FiltPulse[8] ] [#] [0.5 T]£, 

^Select [Table [i, {i, 1, 2*8}] / / Map [ConvertToBitSeq [8] , #]&, (#[[5]] == 1)&] ] 

' l n'li 92S f' °* 74922 ' 0.749173, 0.749219, 0.711482, 0.711529, 0.711482, 0.711435 
0.829796, 0.82975, 0.829703, 0.829749, 0.792012, 0.792059 0 792012 0 7919GS 
M^??!' °' 776839 < 0.776791, 0.776838, 0.739101, 0.739147, V^sI^O 73^4 ' 
0.802178, 0.802131, 0.802084, 0.80213, 0.764394, 0.76444, 0.764393 0 764347 
M^f^' M 59475 ' °* 759428 ' 0.759474, 0.721738, 0.721784, 0.721737, 0 72169 
n*^^; 1 '. 0 ;^ 9495 ' °* 8194 47, 0.819494, 0.781757, 0.781804 0.781756 0 78171 
0.78714, 0.787094, 0.787047, 0.787093, 0.749356, 0.749403, 0.749355 0 749309 
I'llllll' °* 791876 ' °* 79 ^29, 0.791875, 0.754138, 0.754185, 0.754138, 6.754091 
0.749266, 0.74922, 0.749173, 0.749219, 0.711482, 0 . 71152 9 , 0 . 7 11482 , 0 711435 
0.829796, 0.82975, 0.829703, 0.829749, 0.792012, 0.792059, 0 . 792012 / 0 ' 791965 
0.776885, 0.776839, 0.776791, 0.776838, 0.739101, 0.739147, 0.7391, 0 739054 
0.802178, 0.802131, 0.802084, 0.80213, 0.764394, 0.76444, 0.764393 0*764347 
0.759521, 0.759475, 0.759428, 0.759474, 0.721738, 0.721784, 0.721737 0 72169 
0.819541, 0.819495, 0.819447, 0.819494, 0.781757, 0.781804 0.781756 0 78171 
0.78714, 0.787094, 0.787047, 0.787093, 0.749356, 0.749403, 0.749355 0 749309 
0.791922, 0.791876, 0.791829, 0.791875, 0.754138, 0.754185, 0.754138, 0.754091} 

Hap[(l/2 - cr/2Erf [-i]j& 7 xl] // Apply [Plus, #]&//# /128 &; 
Ber[a_] : = Evaluate [%4 8] 
Ber [0.01] 

0.495 

Erf [Infinity] 
1 

D[Erf [x], x] 

2 E-* 2 



? Erf 



Erf [z] gives the error function erf (z) . Erf[z0, zl) gives the generalized error 
function erf(zl) - erf(zO). 

Plot [Erf [x], {x, 0, 10}] 
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0.9999 
0.9998 
0 . 9997 
0.9996 

Graphics 



? FiltPulse 
Global *FiltPulse 
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T := 3 / 812500 
Null 

?? T 

This is the symbol period 
?? T . 
This is the symbol period 
T : = 3 / 812500 

Clear[T, xO, xl, x2, x3, x4, x5, x6, x7] 

j^.i^s.j^aso ) Pulse[0] [ <5T] + j< G «-* + a H-5 * Pulse[0][2T + 6T] + 

CT°«-* Pulse [0] [4 T + <5T] + tf<-°«-s> Pulse [0] [6 T + <5T] + 

j(-a N . 5 -a M . t -o M . 7 ) p u i s ©[0]t8T + <5T] + j£°m-* + ^-3 p u lse[l][T + 5T] + 

j(on- 3 ) Pulse[l][3T + <5T] + J<- a «- «>' Pulse [ 1] ,[5 T + <5T] / . { a M -+ xO, o H .! -> xl, 
a N - 2 -+ x2, a N _ 3 -+ x3, a N _4 -» x4, c N . 5 -* x5, cr K _ e x6 # cc„_ 7 -+ x7} 

(1. DxO*xl4 X 2*x3*x4 p ulse[ 0] [<5T] + (1. I)^ + x3.x4 p ulse[0 j [2 T + c5T] + 

(1. I) x < Pulse[0] [4T + (5T] + (1. I)- x5 'Pulse[0] [6 T + 6T] + (1. i)-xs-x6-x? Pulse[0 ] [8T + 6T] + 
(1. i)*i**3*** Pulse[l] [T + <5T] + (1. I)* 3 Pulsejl] [3 T + <5T] + (1. I)" x6 Pulsefl] [5 T + 6T) 

jCaw-s) p u lse[l] [3 T + <5T] /. { a N -+ xO, a„ _ 2 -» xl, a w . 2 x2, a H _ 3 -> x3, 
a„. 4 -+ x4, a N . 5 -> x5, a N _ 6 -+ x6, a K _ 7 -*x7} 

(1. I) x3 Pulse [1] [3T + aT] 

(1. » I) x3 Pulse[l] [3 T + <5T] 
(1. I) x3 Pulse[l] [3T + <5T] 
? a 




J = 



